Ghadimi et al. BMC Gastroenterology 201 1, 1 1:104 
http://www.biomedcentral.eom/1 471 -230X/1 1 /1 04 



Gastroenterology 



RESEARCH ARTICLE Open Access 



Family history of the cancer on the survival of 
the patients with gastrointestinal cancer in 
northern Iran, using frailty models 

Mahmoodreza Ghadimi, Mahmood Mahmoodi" Kazem Mohammad, Hojjat Zeraati, Mahboobeh Rasouli and 
Mahmood Sheikhfathollahi 



Abstract 

Background: Gastrointestinal (Gl) tract cancer is one of the common causes of the mortality due to cancer in 
most developing countries such as Iran. The digestive tract is the major organ involved in the cancer. The northern 
part of the country, surrounded the Caspian Sea coast, is well known and the region with highest regional 
incidence of the Gl tract cancer. In this paper our aim is to study the most common risk factors affecting the 
survival of the patients suffering from Gl tract cancer using parametric models with frailty. 

Methods: This research was a prospective study. Information of 484 cases with Gl cancer was collected from Babol 
Cancer Registration Center during 1990-1991. The risk factors we studied are age, sex, family history of cancer, 
marital status, smoking status, occupation, race, medication status, education, residence (urban, rural), type of 
cancer, migration status (indigenous, non-native). The studied cases were followed up until 2006 for 15 years. 
Hazard ratio was used to interpret the death risk. The effect of the factors in the study on the patients survival are 
studied under a family of parametric models including Weibull, Exponential, Log-normal, and the Log-logistic 
model. The models are fitted using with and without frailty. The Akaike information criterion (AlC) was considered 
to compare between competing models. 

Results: Out of 484 patients in the study, 321 (66.3%) were males and 163 (33.7%) were females. The average age 
of the patient at the time of the diagnosis was 59 yr and 55 yr for the males and females respectively. 
Furthermore, 359 (74.2%) patients suffered from esophageal, 110 (22.7%) patients recognized with gastric, and 15 
(3.1%) patients with colon cancer. Survival rates after 1, 3, and 5 years of the diagnosis were 24%, 16%, and 15%, 
respectively. We found that the family history of the cancer is a significant factor on the death risk under all 
statistical models in the study. The comparison of AlC using the Cox and parametric models showed that the 
overall fitting was improved under parametric models (with and without frailty). Among parametric models, we 
found better performance for the log-logistic model with gamma frailty than the others. Using this model, gender 
and the family history of the cancer were found as significant predictors. 

Conclusions: Results suggested that the early preventative care for patients with family history of the cancer may 
decrease the risk of the death in the patients with Gl cancer. The gender appeared to be an important factor as 
well so that men experiencing lower risk of death than the women in the study. Since the proportionality 
assumption of the Cox model was not held (p = 0.0014), the Cox regression model was not an appropriate choice 
for analysing our data. 
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Background 

Cancer is known as one of the major causes leading to 
many disorders, death, and disabiUties worldwide [1,2]. 
Cancer has affected increasingly the human population 
during the past decades so that considerable amount of 
health care resources have been allocated to diminish its 
side effects [3]. It is predicted to become the leading 
cause of death in many developed and developing coun- 
tries such as Iran [1,4]. 

Esophageal, Stomach, and Colorectal cancers are three 
most common types of cancer among Iranian popula- 
tion [5]. The Northern part, located on the Caspian 
coast has been reported as the main area of the county 
dominated by Gastrointestinal (GI) Malignancies tract 
cancer [6,7]. Out of 70.4 million of the country's popula- 
tion, almost 50,000 new cases of cancer are reported 
each year. In more than 38% cases, the GI tract is par- 
tially or completely affected. 

Stomach, esophageal, and colon cancer are three most 
common types of cancers reported among males. For 
females, the breast cancer should be added to this list 
[8,9]. Cancer is the third most common cause of death 
in Iran. This accounts for 14% of the total mortality of 
which the GI cancer accounts for approximately half 
(44.4%) of all cancers related deaths [10,11]. Unfortu- 
nately, the GI cancer in Iran is diagnosed when the dis- 
ease is in its developed phase and hence the patience 
the available therapies treatment are less effective to 
cure the patience [1,11]. Practically, the early diagnosis 
of the GI tract cancer gives more chance to the patients 
to recover from the discomfort. 

Survival data are often modeled using the Cox propor- 
tional hazards model which estimates the covariate 
effects as the log hazard ratios. This model is free of 
estimating the baseline hazards for the model. However, 
since the hazard function is directly related to the time 
course of the disease, its behavior may be of medical 
interest. The baseline hazard rate can help us to under- 
stand the common history of the disease respect to the 
hazard rate changing over time [12,13]. Coxs semi-para- 
metric regression model [14] is frequently used to ana- 
lyse the survival data. Alternatively the fully parametric 
models such as Weibull, Log-Logistic and Log-Normal 
models can be used [15,16]. They can offer a gain that 
may not be obtained under Cox s model. Efron [17] and 
Oakes [18] showed asymptotically that under certain cir- 
cumstances, parametric models can lead to more effi- 
cient estimates of the parameter. 

In survival analysis, to model the data in which the 
mortality reaches a peak and then starts to decline, a 
model with a non-monotonic (hump-shaped) failure rate 
can be used. This is the case with our data we use in 
this paper. In order to capture efficiently this property 



of data, the Log-logistic and Log-normal model are 
often used [13]. However, If the issue of outliers is not 
major, the Log-logistic model can be used to approxi- 
mate the Log-normal model. Moreover, for censored 
data, the Log-logistic model has a simple hazard form 
and survival function [12,13]. For these reasons, we use 
the Log-logistic function in this paper for analyzing our 
data. 

The aforementioned pattern for hazard function was 
the case in our study. Hazard function increased slowly 
until after a while started to decline. Because of this pat- 
tern in our data Cox, Weibull, and Exponential models 
are not appropriate ones and as was said in above Log 
logistic model seems better as results of our findings 
verified the issue. 

It is assumed that for the unique covariate inputs, the 
survival function under the Cox proportional hazards 
and parametric models is the same for subjects. How- 
ever, the data may present extra-variation due to the 
unobserved factors. In this study, we collected data on 
all possible factors we thought might influence the 
patient's survival. 

A model becoming increasingly popular for modeling 
the multi-level individual survival times is frailty model. 
A frailty is an unobserved random effect shared by sub- 
jects within a sub-group. Frailty models are also used to 
capture the overdispersion in univariate survival studies. 
In this paper, the frailty refers to the effect of the unob- 
served factors on the subject's survival. Ignoring frailty 
may lead to the biased survival estimates. The overdis- 
persion is modeled using a latent multiplicative effect 
on the hazard, or frailty. A gamma or inverse-Gaussian 
distribution is commonly used to model the frailty 
[12,19-22]. Thus, the hazard of a population is inter- 
preted as the mean of individual hazards among the sur- 
vivors. Frail individual with notable values of frailty will 
tend to die sooner [19]. The frailty (random effect) can 
be integrated out (in closed form or by numerical or 
stochastically integration, depending on the frailty distri- 
bution) to get a likelihood function not depending on 
unobserved quantities [19]. 

By the expectation is conditional on being at risk at 
time point t, it mention averaging over a subset of the 
original population. Therefore, relative weights for 
hazards with high frailty become smaller as time goes 
by, corresponding to high mortality. An important 
implication is that studies of human aging based on 
cohort mortality data may be systematically biased or 
based on erroneous functional forms [19]. The aim of 
the this paper is to investigate the factors influencing 
the survival of the patients with GI tract cancer using 
parametric models with frailty. We also compare our 
results with that of achieved under the Cox model. 
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Methods 

This survey was a prospective study. The total number 
of 484 patients with developed GI tract cancer registered 
at the Babol Cancer Registration Center during 1990- 
1991. They then followed up for 15 years until 2006. 
The socio-demographic and clinical data obtained using 
questionnaire and the patients' clinical records. Written 
informed consent from patients was obtained prior to 
entering the study. Patients completed a questionnaire 
that assessed satisfaction with the informed consent pro- 
cedure. Also to maintain patient privacy, all records 
were coded with a unique project identifier prior to 
transmission to the data collection. The study was con- 
firmed by the Ethics Committee of Tehran University of 
Medical Sciences. The factors we consider in our study 
are age at diagnosis, gender, place of residence, pro- 
vince, type of cancer, method of cancer detection, family 
history of cancer, education, job, marital status, cigarette 
smoking, ethnicity, migration status, drug use. 

A multivariate parametric regression model (with and 
without frailty) was developed to analyse the prognostic 
factors related to the longevity of patients. To compare 
the different parametric models and their efficiency the 
Akaike Information Criterion (AIC) [23], Cox-Snell, and 
deviance residual plots were used. The AIC was consid- 
ered to assess the general goodness of fit of the statisti- 
cal models. The lower value of the AIC, the better 
model to fit the data. Hazard rate (HR) [12,24] was used 
to interpret the death risk of the parametric models. For 
the statistical analysis, the statistical software SAS 9.1 
and STATA 8.0 were used. The values less than 0.05 for 
probability, p < 0.05, was defined as the level of our sta- 
tistical significance. 

Results 

Out of 484 initial patients with developed GI cancer, 
321 (66.3%) were men and 163 (33.7%) women. The 
mean ± standard deviation of age at diagnosis was 58.26 
± 10.9 years and the median survival time was found 9.1 
months. The estimated survival rates in 1, 3, and 5 years 
after diagnosis were 0.24, 0.16, and 0.15 respectively. 
The type of cancer in these patients was as follows: eso- 
phageal (74.2%), stomach (22.7%) and colon (3.1%) 
(Table 1). During the following up, the total number of 
426 (88.0%) deaths were observed (non-censored obser- 
vations) and 88 (12.0%) patients survived or exact details 
of their survival status were not available (Loss to follow 
up) (right censored observations). 

According to the fact that the proportionality assump- 
tion of Cox model was not met in our data (p = 
0.0014), using Cox regression was not suitable, even 
adding frailty term (with gamma and inverse-Gaussian) 
in to Cox model, proportionality assumption was ever 
violated and there was no remedy in the violation of the 



Table 1 Characteristics of patients with Gastrointestinal 



tarct cancer diagnosis 


Characteristic 


n (%) 


Gender 




Male 


321 (66.3) 


Female 


163 (33.7) 


Place of residence 




Rural 


256 (52.9) 


Urban 


228 (47.1) 


Province 




Mazandaran 


288 (59.5) 


Golestan 


196 (40.5) 


Type of cancer 




Esophageal 


359 (74.2) 


Stomach 


110 (22.7) 


Colorectal 


15 (3.1) 


Method of cancer detection 




Clinical diagnosis 


35 (7.2) 


Direct endoscopy and biopsy 


39 (8.1) 


Conventional chest x-ray 


410 (84.7) 


Family history of cancer 


142 (29.3) 


Education 




Literate 


52 (10.7) 


Illiterate 


432 (89.3) 


Job 




Farmer 


252 (52.1) 


Employee 


7 (1.4) 


Others 


225 (46.5) 


Marital status 




Married 


459 (94.8) 


Single 


25 (5.2) 


Cigarette smoking 


215 (44.4) 


Ethnicity 




Aryan 


327 (67.6) 


Gilak 


12 (2.5) 


Torkaman 


100 (20.7) 


Others 


45 (9.2) 


Migration status 




Native 


430 (88.8) 


Non-native 


54 (11.2) 


Drug use 


119 (24.6) 



PH assumption. Thus Cox model was omitted from 
study. 

The Kaplan-Meier estimates of the survival functions 
for the gender and the family history of the cancer are 
given in the Figure 1. 

Figures 2, 3 plots the Cox-Snell and deviance residuals 
under the parametric models; log-normal, log-logistic, 
and WeibuU model. In overall, the plots show smaller 
residuals using parametric models and therefore we may 
conclude they have better performance than the Cox 
model. Furthermore, the parsimonious of the Cox-Snell 
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2000 4000 
analysis time 



(c) 

Figure 1 Survival curve of Gl tract cancer patients using Kaplan-Meier method, (a), (b) Kaplan-Meier estimates of tine survival curves for Gl 
tract cancer data separated by family history of cancer and gender, respectively, (c) Kaplan-Meier overall survival curves. 



residuals under the log-normal and log-logistic model 
with gamma frailty to the 45 degrees line in Figure 3 
confirms these models provide better fitting to our data. 
It can be also seen that the log-logistic model has better 
performance over the log-normal model. The weak per- 
formance of the WeibuU model which assumes the pro- 
portional hazards can be due to the violation 
assumption of the proportional hazards. 

The similar conclusion can be obtained by using AIC. 
The AIC of each model in the study is given in Table 2. 
The best scores are achieved under the log-logistic 
model. The Weibull model is the next best model fol- 
lowed by the log-normal. Table 2 also suggests the log- 
logistic with gamma frailty as the most efficient model 
for our data. 

Table 3 reports the detailed results of the multivariate 
analysis for the parametric models with and without 
frailty based on the HR for each variable. Results of the 
multivariate analysis show that the family history of the 
cancer appears a significant factor in all fitted models. 
This implies that patients with the family history of the 
cancer are less survived than others. 



Gender is significant under the log-normal and log- 
logistic with gamma frailty model but not significant 
factor under other models. This indicates that the level 
of the death risk due to Gl cancer was reduced signifi- 
cantly for the women in the study during the following 
up period. 

None of the parametric models suggests age, resi- 
dence, province, type of cancer, methods of cancer diag- 
nosis, educational level, occupation, smoking, ethnicity, 
migration status and drug use as a significant prognostic 
factors. 

Discussion 

Gl tract cancer is one of the most common types of 
cancer in Iran [10]. The cancer is a particularly devastat- 
ing form of cancer with a relatively low survival rate, 
and people generally will not live a long time after diag- 
nosis. Several factors known in various studies as influ- 
encing prognosis factors and have been introduced 
[25-33]. 

In the literature, there are many studies on the field of 
cancer, but researchers tend to examine the effects of 
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Log-logistic Model 



Log-Normal Model 



Mean Deviace Residual = 0.184 



V. 

•• • X* % 



Mean Deviace Residual = 0.238 



Figure 2 Deviance residuals to evaluate model fit of parametric models. In this plot, the deviance residual is large for short survival times 
and then decreases with time. This pattern suggests that the log normal and log logistic models with gamma frailty are better than other both 
models (The log logistic model has the lowest mean deviance residual with respect to other models). 



covariates on patients survival using Cox regression 
model instead of parametric ones. A systematic study on 
Cancer Journals shows that only in 5% of studies of can- 
cer in which Cox regression model is used the assump- 
tions of the model have been investigated [34]. If 
presumptions are not met, results of Cox model are ser- 
iously under question. As an alternative, parametric 
models such as log-normal, log logistic, WeibuU, and 
exponential can be employed. The only assumption of 
parametric models is that the variable time follows a 
specific distribution [13,24]. 

In this paper we aimed to study the possible relation- 
ship between the survival of the patients with GI tract 
cancer and several most common prognosis factors such 
as age at diagnosis, gender, place of residence, province, 
type of cancer, method of cancer detection, family his- 
tory of cancer, education, job, marital status, cigarette 
smoking, ethnicity, migration status, drug use. 

We found gender and the family history of the cancer 
significant prognostic factor. This supports the past stu- 
dies reporting better survival for women with developed 
GI tract cancer [1,35-37] and the family history of the 
disease as a significant factor [38,39]. 



Statistical assessment using AIC of the studied models 
showed that the log-logistic model with gamma frailty 
will describe our data better. 

Due to their better performance, our intension is to 
use parametric models. However, the efficiency of the 
parametric models is greatly affected by the volume of 
the censored observations. For having sensible results, it 
is recommended that the percentage of the censored 
data should not be more than 50% [40]. This condition 
is satisfied with our data as they consist of 15% censored 
observations. 

Nardi et al. [40] compared the performance of the 
Cox model and some parametric models. They used 
normal-deviate residuals [41] to evaluate the assump- 
tions of parametric models. They also studied WeibuU 
model based on the estimated variation of the parameter 
rate criteria, and concluded that the Weibull was the 
superior model. In our study, we found the log-logistic 
model to have better performance than the other mod- 
els in the study. 

By a simulation study, HRbe et al. [42] compared the 
Cox regression model and the accelerated failure time 
(AFT) models. They used the proposed method by Stute 
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Residuals 



Weibull 



Residuals 
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Residuals 



Residuals 



Figure 3 Cox-Snell residuals obtained from fitting various survival models with gamma frailty to the Gl tract cancer data. The panels 
indicate the Cox-Snell residuals (together with their cumulative hazard function) obtained from fitting different parametric models to the same 
data via maximum likelihood estimation. 



Table 2 Overview of the Akaike Information Criterion 
Scores 





AlC 


Rank 


Without Heterogeneity 






Exponential 


2127.24 


4 


Weibull 


1918.18 


3 


Log-normal 


1800.44 


2 


Log-logistic 


1786.94 


1 




AlC 


Rank 


Gamma Heterogeneity 






Exponential 


1780.94 


4 


Weibull 


1739.86 


2 


Log-normal 


1746.16 


3 


Log-logistic 


1717.84 


1 




AlC 


Rank 


Inverse Gaussian Heterogeneity 






Exponential 


1851.98 


4 


Weibull 


1851.58 


2 


Log-normal 


1759.78 


3 


Log-logistic 


1725.02 


1 



[43] for fitting linear regression models with right-cen- 
sored data. Their results showed that whatever the pro- 
portional hazards assumption is violated or not, the log 
logistic, log-normal, and the Stute models are more effi- 
cient than the Cox model. 

Bradburn et al. [44] evaluated the adequacy of some 
parametric models and the Cox proportional hazards 
model using model's residuals and the AIC. They found 
that the generalized gamma model and parametric mod- 
els achieved both a higher log-likelihood and a lower 
AIC. 

For the Cox and parametric models, the hazard func- 
tion may depend on the unknown or latent factors 
which can lead to the biased estimates of the regression 
coefficients [19,45]. To overcome this issue we used the 
frailty models. In fact these models are used to explain 
the random variation of the survival function that may 
exist due to unknown risk factors such as genetic factors 
and other environmental factors [19,22,45-47]. 

Random effects models are known as the frailty mod- 
els in the survival analysis. These models, widely studied 
in the 1990's, are relatively new in the survival field and 
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Table 3 Multivariate model of parametric models with and without frailty 



vdridbl6S 




Expon6ntidl 






Weibull 

up (p i/aliio^ 






LoQ-normal 

up m i/aliio^ 






LoQ-loQistic 






without 


Gamma 


Inv- 


without 


Gamma 


Inv- 


without 


Gamma 


Inv- 


without 


Gamma 


Inv- 




frailty 


frailty 


Gussian 


frailty 


frailty 


Gussian 


frailty 


frailty 


Gussian 


frailty 


frailty 


Gussian 








frailty 






frailty 






frailty 






frailty 


Age(years) 


1.005 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 


1.01 




[{j.ZyO) 


i^u.zy 1 ) 


l,U.z44j 


(r\ ^1 


[U. \ / Z) 


XSJ.Z^i) 


{{j.ZjZ) 


yj.Zoo) 


ir\ lA Q^ 

iU.Z4oj 


[u.zoy) 


iU.zuoj 


[u. \ oo) 


Gender 


1.12 


1.32 


1.25 


1 .1 1 


1 .47 


1 .22 


1.37 


1 .64 


1.61 


1.28 


1.45 


1 .47 






[[J.Z5Z) 


l,U.oz4j 


iU.4Dyj 


iu.ujyj 


IC\ Q/l 


[{JA^Z) 


/'nni 1 ^-^ 
iU.U 1 1 1 


[{J.UZZ) 


/n lA 
iU.Z4jj 


(r\ n/i 1 
l,U.U4 1 ) 


(r\ n/i 
l,U.U4Zj 


Place of 


1.15 


0.99 


1 .07 


1 .07 


0.95 


1.06 


1 .02 


0.96 


0.98 


0.99 


0.92 


0.93 


residence 






l^U.DDoj 


(Pt c:^ A\ 
[{J.J 1 4j 


iU./ I4j 


iu.oy/j 


(r\ Qni ^ 
iu.yu 1 1 


yj./Oo) 


fr\ QQ1 ^ 

iu.oy 1 ) 


yj.yyZ) 


[UAyZ) 


l^U.JDOj 


Province 


0.76 


0.96 


0.77 


0.79 


1 .09 


0.79 


0.87 


1 .06 


1 .01 


0.93 


1 .12 


1.12 




[U.UOy) 


(0 9.1 A\ 


\\J. \ 5Z) 


[U.UDU) 


VU.D 1 1 ) 


\y).\o\) 


(C\ A1 1\ 
l,U.4 1 Z) 


yj./yjy) 


(0 QQQ\ 
[U.yOy) 


yJ.OOy) 


[U.DyD) 


VU.4JOj 


Type of cancer 


0.87 


0.95 


0.88 


0.9 


0.99 


0.88 


0.93 


1 .05 


1 .02 


0.94 


1 .03 


1 .02 




[u.ZUy) 


[[j./y/) 


(n A 

l,U.4zUj 


[{j.dZo) 


iu.yooj 


l,U.43Zj 


yj.oZj) 


yj./ Zo) 


iu.yz 1 ) 


^'n 71 


[{J.ODJ) 


(c^ Qni ^ 
l^u.yu 1 ) 


Family history of 


1 .85(< 


1 .69 


1.92 (< 


1 .54 (< 


1 .49 


1.82 (< 


1 .82 


1 .47 


1 .56 


1 .72 


1 .39 


1 .43 


cancer 


U.UU 1 ) 


iU.UU4j 


U.UU 1 ) 


U.UU 1 ) 


l^u.uuyj 


n nni \* 
U.UU 1 ) 


ic\ nni ^-^ 
1,U.UU 1 ) 




l,U.UU/j 


(n nni ^^ 
iU.UU 1 ) 


[u.u \ o) 


iU.U 1 D) 


Education 


0.6 


0.63 


0.6 


0.69 


0.74 


0.62 


0.62 


0.73 


0.7 


0.6 


0.79 


0.77 






[U. \ \ 4j 


iU.Uj4j 


iU.UoOj 


\SJ.Z \ 5) 


^'n nc^Q^ 
l^u.Ujyj 


iU.U/ 1 ) 


iU.ZU 1 ) 


/'n 1 c^n^ 


(n nAn^ 

l^U.UDUj 


[[J.DU/ ) 


(n new 
iU.Z/Uj 


Marital status 


1.61 


1.64 


1.75 


1.45 


1.49 


1.69 


1.67 


1.35 


1.47 


1.67 


1.31 


1.39 




(0.063) 


(0.224) 


(0.139) 


(0.156) 


(0.222) 


(0.150) 


(0.180) 


(0.379) 


(0.272) 


(0.177) 


(0.361) 


(0.305) 


Cigarette 


1.35 


0.99 


1.23 


1.2 


0.9 


1.2 


1.06 


0.84 


0.9 


1.03 


0.85 


0.87 


smoking 


(0.014)^ 


(0.989) 


(0.263) 


(0.109) 


(0.533) 


(0.289) 


(0.718) 


(0.306) 


(0.528) 


(0.871) 


(0.298) 


(0.363) 


Migration status 


1.40 


1.35 


1.43 


1.23 


1.37 


1.39 


1.45 


1.35 


1.41 


1.35 


1.28 


1.3 




(0.047)^ 


(0.266) 


(0.154) 


(0.194) 


(0.194) 


(0.170) 


(0.157) 


(0.191) 


(0.157) 


(0.227) 


(0.267) 


(0.238) 


Drug use 


0.95 


1.11 


1.01 


0.97 


1.07 


1.01 


1.06 


1.11 


1.1 


1.11 


1.09 


1.09 




(0.727) 


(0.624) 


(0.958) 


(0.837) 


(0.674) 


(0.944) 


(0.758) 


(0.561) 


(0.614) 


(0.598) 


(0.618) 


(0.637) 


Method of cancer 


0.64 


0.89 


0.73 


0.76 


1.11 


0.76 


0.84 


1.11 


1.06 


0.85 


1.07 


1.08 


detection (Direct 


(0.023)* 


(0.751) 


(0.292) 


(0.259) 


(0.814) 


(0.319) 


(0.573) 


(0.712) 


(0.839) 


(0.613) 


(0.790) 


(0.775) 


endoscopy and 


























biopsy) 


























Method of cancer 


0.39 (< 


0.57 


0.44 


0.51 


0.70 


0.46 


0.56 


0.95 


0.84 


0.52 


0.93 


0.88 


detection (Clinical 


0.001)* 


(0.142) 


(0.01 1)* 


(0.003)* 


(0.495) 


(0.013)* 


(0.071) 


(0.878) 


(0.586) 


(0.059) 


(0.795) 


(0.688) 


diagnosis) 


























juu vi^i 1 1 piuy trey 


1 06 


0 85 


0 95 


0 98 


0 69 


0 93 


1 1 3 


1 89 


1 71 


0 85 


0 84 


0 83 




(0.895) 


(0.829) 


(0.937) 


(0.972) 


(0.758) 


(0.917) 


(0.856) 


(0.445) 


(0.494) 


(0.814) 


(0.773) 


(0.772) 


Job (Others) 


1.01 


0.98 


1.01 


0.99 


0.99 


1.01 


1.04 


1.08 


1.09 


0.97 


0.99 


0.99 




(0.926) 


(0.918) 


(0.944) 


(0.996) 


(1.00) 


(0.972) 


(0.832) 


(0.638) 


(0.618) 


(0.873) 


(0.936) 


(0.967) 


Ethnicity (Gilak) 


1.35 


0.53 


0.87 


1.08 


0.22 


0.85 


0.62 


0.39 


0.42 


0.58 


0.43 


0.42 




(0.329) 


(0.211) 


(0.783) 


(0.806) 


(0.041)* 


(0.745) 


(0.351) 


(0.027)* 


(0.061) 


(0.224) 


(0.020)* 


(0.023)* 


Ethnicity 


1.33 


1.30 


1.35 


1.23 


1.34 


1.32 


1.24 


1.04 


1.04 


1.34 


1.23 


1.21 


(Torkaman) 


(0.089) 


(0.327) 


(0.255) 


(0.202) 


(0.434) 


(0.255) 


(0.414) 


(0.878) 


(0.881) 


(0.231) 


(0.293) 


(0.351) 


Ethnicity (Others) 


1.18 


1.39 


1.21 


1.12 


1.79 


1.21 


1.23 


1.26 


1.21 


1.39 


1.52 


1.48 




(0.377) 


(0.267) 


(0.490) 


(0.517) 


(0.160) 


(0.476) 


(0.469) 


(0.344) 


(0.465) 


(0.229) 


(0.062) 


(0.092) 



* Significant at 0.05 level 
HR, Hazard rate 



are currently mainly under investigations, but technical 
problems in estimating the parameters of frailty models 
made to be used less compared to the Cox model. 
Using frailty to model the extra-variation in univariate 
lifetime data goes back to the work of Vaupel et al. 
[48]. Henderson and Oman [49] in a theoretical method 
revealed that in case of non-use of frailty model when 
there is frailty effect bias may occur in the estimates of 
regression coefficients. Schumacher et al. [50] showed 
that ignoring an important factor can lead to lower 



estimations of the relative risk by the fitted models. 
Keiding et al. [51] showed how removing one of the 
two explanatory variables might increase the variance of 
the hazard function and biased estimation of other coef- 
ficients in the fitted model. They suggested using AFT 
models to handle the effect of unobserved variables. 

According to our findings, log logistic model with 
gamma frailty is more suitable statistical model in survi- 
val analysis in patients with GI cancers rather than 
other parametric models. 
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Conclusions 

Our study showed that the gender and the family his- 
tory of the cancer were two factors that can significantly 
affect the lifetime of the patients with GI tract cancer. 

According to our findings the early recognition of 
family history of cancer and, in consequence, awareness 
of family members to consider the possibility of family 
screening may result in a decrease in death rate due to 
GI tract cancer. 

Furthermore, we found that the death risk of the GI 
tract cancer for the men was significantly lower than the 
women. We also recommended to use the log-logistic 
with gamma frailty model, to evaluate the effects of the 
prognostic factors on the developing the GI tract cancer. 

Limitation 

One of the limitations of this study was the lack of an 
efficient recording medical system in the Babool Cancer 
Registeration Center. Currently there is no any informa- 
tion available for some clinical factors such as the type 
of esophageal cancer (adenocarcinoma, squamous) and 
the stage of the disease. 
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